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OO Abstract: Any careful singularity analysis of the Mixmaster Universe uncovers the instance 
^ of a triple —1 resonance. The Mixmaster Universe does not exhibit a closed-form solution 
CN . and so a correct interpretation of the meaning of the triple —1 resonance is difficult. We 
provide a system of differential equations for which both a closed-form solution and a triple 
— 1 resonance exist. 
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1 Introduction 



^ ! Our understanding of the workings of a physical system is generally enhanced by mathe- 
\Q \ matical analysis provided that analysis is undertaken in a state of confidence in the meaning 
^ ; of the results obtained. One way to develop such confidence is to have experience of the ap- 
• ■ plication of that analysis to a model system of known properties. This is not a perfect route 
^ ■ since in exploring the unknown one may stray so far from the known that prior experience 
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00 , cannot help. In this Letter we intend to highlight this problem in terms of the singularity 
analysis of ordinary differential equations. In Cosmology the Bianchi Type IX model, com- 
_ ^ monly known as the Mixmaster Universe, has received a number of investigations through 
^ singularity analysi^ O [71 [H [U [TOl [151 M [HI [2] • The conclusions reached reveal an interesting 
evolution in the various interpretations of the results of the analysis. The first paper [9j^l 
reported a four-parameter particular solution. Initially Contopoulos et al concluded that 
the system of differential equations describing the Mixmaster Universe satisfied the Painleve 
criterion for integrability. Subsequently Contopoulos et al moved to the position that the 
equations were not integrable due to the presence of essential singularities. This opinion has 
been reinforced by the conclusions of Demaret and Scheen, Latifi et al and Springael et al. 
More recently the present authors proposed that there existed two patterns of singularity be- 
haviour, if one omits cyclic permutations, both of which possess the full number of constants 



^permanent address: School of Mathematical Sciences, University of KwaZulu-Natal, Private Bag X54001 
Durban 4000, Republic of South Africa 

^There have been other approaches and the reader is referred to [2] for a recent summary of the literature. 

•^Although the first of the three papers by Contopoulos and his coworkers appeared earlier in the Journal 
of Physics A: Mathematical and General, the paper by Cotsakis et al was submitted over a month earlier to 
the same journal, but it did not appear until some seven or eight months later. 
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of integratior0 [2]. 



It is certainly not unreasonable to wonder why various investigations of the same system 
have led to a variety of interpretations. There appear to be three sources for concern. 

The first is the understanding of the meaning of negative resonances - apart from the 
generic —1 - arising from the analysis. In the case of an equation possessing only nonpositive 
resonances Lemmer et al [17] had already demonstrated that one could have a Laurent 
expansion with the powers commencing at the singularity and descending to minus infinity. 
This was formalised by Feix et al [I3]. The interpretation of such a series is that it is valid 
on the exterior of a disc centred on the singularity so that it has the nature of an asymptotic 
expansion. 

The second is concerned with the existence of both positive and negative resonances - 
again apart from the generic —1 - for a given leading-order behaviour. Latifi et al proposed 
that these represented two particular solutions, one descending and one ascending, since the 
number of constants of integration for each series would be less than the number required 
|T5]. In a recent paper [1] it was demonstrated by way of explicit example that the case of 
resonances of mixed sign was simply a Laurent series valid on an annulus centred on the 
singularity. 

The third source for concern is the appearance of a triple —1 resonance for one of the 
branches of the equations of the Mixmaster Universe. Originally Cotsakis et al [9j thought 
that it indicated a four-parameter particular solution. Then Contopoulos et al |6] proffered 
the opinion that it was due to the forcing of the three functions to have a singularity at 
the same point. In our more careful analysis - having benefited from the evolution of ideas 
over the past decade - we found that the triple —1 resonance admitted three constants of 
integration [2]. Indeed the Laurent series contained the full set of constants of integration 
and would exist on an annulus since there were both positive and negative resonances. 
There were two Laurent expansions. One was moderately conventional and would exist on a 
punctured disc centred on the singularity and then there was the one on the annulus. This 
result melds very nicely with the numerical studies of Bountis and Drossos [S]. A matter 
of particular interest in this case is that the presence of the triple —1 resonance indicates a 
Laurent series of this type. It is a commonplace that the existence of a resonance at —1 is 
generic and is related to the location of the movable singularity. That was contradicting the 
whole meaning of the resonances which were supposed to indicate the first entrance of yet 
another arbitrary constant in the expansion pTj. The reasoning strengthens when we move 
to systems of differential equations and that a multiple —1 resonance can indicate additional 
constants of integration whereas for a single equation it is not good news. 

In this Letter it is our intention to discuss the singularity analysis of a system of three 
second-order nonlinear ordinary differential equations of moderately complex form which 
nevertheless can be integrated explicitly. Thus we can see how the various principal branches 

"^Since the differential equations of the model are three of the second order, this means six. As there is a 
constraint, the number is reduced to five. Since the constraint is a particular value of a first integral of the 
system of differential equations, it plays no role in the analysis except to reduce the solution to existing on 
a surface in the space of initial conditions. 
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relate to the general solution. 



2 The System and Its Solution 

We consider the system 



X + 3 {xx + yz + zy) + + + + 6xyz = 

y + 3 {xy + yx + zz) + 3 (^x'^y + y'^z + z'^x^ = (2.1^ 

^ + 3 {xz + yy + zx) + 3 ixy"^ + yz^ + zx^ ] = 



which has the three Lie point symmetries 

Ti = dt, T2 = -tdt + xd^ + ydy + zd, and = -t^dt + 2{-l + tx)d^ + 2tydy + 2tzd, (2.2) 
with the algebra si {2, R). 

The solution of the system (12. ip is easily found to be 

C^A + A^B + B'^C - ABA - BCB - ACC 
^^^^ " A3 + S3 + c'3 - 3ABC ^^-^^ 

B'^A + C'^13 + A^C - AC A - ABB - BCC 
~ A3 + S3 + C3 - 3ABC ' 

where the functions A{t), B{t) and C{t) are given by 

A{t) = ao + bot + CQt^ 

Bit) = ai + bit + cit^ (2.4) 
C{t) = a2 + M + C2t^ 

and are introduced to render the expressions in (12.31) to be of a manageable size. Although 
there are nine arbitrary constants listed in (12. 4p . only six in the solutions, (12.30 . are in fact 
independent. 

It is evident from (\2.3\\ and (12.40 that the solution of the system is analytic away from 
its polelike singularities. Thus we can proceed with the singularity analysis in the sure 
knowledge that we can interpret the results in light of our possession of the explicit solution. 
Such explicitly integrable systems may not give a complete answer to every question, but 
they do help to prevent the making of unjustified inferences. 



3 Singularity Analysis 
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Since system f|2.ip possesses the self-similar symmetry, all terms are dominant^! and 
it is obvious from the symmetry that the exponent of the leading-order terms is —1. We 
make the usual substitution, x = aT~^, y = /3t~^ and z = jt~^, where r = t — to, into the 
equations. The resulting set of algebraic equations is 

2a - Sa^ - 6/?7 + + + 7^ + 6a/37 = 

2/3 - 6a/? - 37^ + Sa^Z? + 3/3^7 + 3a72 = (3.1) 
27 - 6a7 - 3/3^ + 3a/32 + 3/37^ + 3a27 = 

which provides a large number of solutions - indeed 27. For each one of these solutions which 
does not contain a zero value for the coefficient of a leading-order term - that would make 
for a particular solution - we can perform the calculations for the resonances. We make the 
substitutions, x = ar'^ + ^t'^'^, y = (3t^^ + ut'^^^ and z = jt^^ + ^t'^~^, where s is the 
value of the resonance at which arbitrary constants enter the Laurent expansion, into the 
full equations since all terms are dominant for this leading-order exponent. We summarise 
our results in Table I. 



TABLE I BE HERE 

Table I: Coefficients of the leading-order terms and corresponding resonances for system 

(12.11) in the case that all terms are dominant. 



The order of presentation of the different possible patterns of values of resonances in 
Table I is dictated by the number of resonances of value —1 which happens to match an 
ordering of the values of a from zero to two by increments of one third. Within each triplet 
we note that the sums of the values of /? or 7 in the second and third members of the triplet 
are the negatives of the value of the same for the first member of the triplet. Doubtless there 
could be other orderings based upon some other criteria. For example one could take the sum 
of the values of a, (3 and 7 and use this to order the sets of the values of the resonances. It is 
quite evident that this system of second-order equations displays an unexpected richness in 
both the coefficients of the leading-order terms and the values of the resonances. In the case 
of certain sequences of scalar equations, such as the Riccati and Emden-Fowler sequences, 
this richness of property has been observed [HI [3l [16] and, if one could construct a sequence 
commencing with the present system by finding its recursion operator, the patterns could 
become even more exotic. 



In our discussion we ignore the sets of coefficients of leading-order behaviour which con- 
tain a zero or more. 

The first triplet gives a standard representation of a Right Painleve Series in that the 
nongeneric resonances are positiv^. 

The second triplet has one of the +1 resonances replaced by a —2. We now have a case 
of resonances of mixed sigE0. 

^Here we do not make an investigation of subdominant behaviours since the analysis of the aU-terms- 
dominant case provides more than sufficient results for our present purposes. 

^The —1 resonance is always to be expected in this analysis [14 . The interpretation of its meaning is 
another matter and is addressed below with this system providing a specific medium. See also 

^Since the —1 resonance is expected, curiously it is not normally taken to be counted in a discussion of 
the existence of resonances of mixed sign. 
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For the third triplet one of the +2 resonances of the first triplet has moved down to — 1 
so that we have a double resonance at —1. We return to this case below. 

The following two triplets are related to the third in the sense that one of the +1 reso- 
nances has become —2. 

In the instance of the sixth triplet a second of the +1 resonances has moved down to the 
value —2. 

At the seventh triplet the remaining +2 resonance has made the third of the —1 reso- 
nances. 

When we come to the eighth triplet, the pattern of resonances is the direct opposite of 
that for the first triplet. 

Corresponding to each value of the resonances for each of the sets of {a,/?, 7} there is 
an eigenvector. To present them for all the possible coefficients of the leading order terms 
would be pleonastic. However, we believe that it is informative to consider the third and the 
eighth triplets. In the case of the third triplet we have the intrusion of a double resonance 
at the so-called generic value —1. For the (seventh and) eighth triplet we have a triple — 1 
resonance. 



For the first member of the third triplet the eigenvectors are (in the order corresponding 
to the ordering of the values of the resonances) 
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where u and ^ are parameters, and for the second 
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The eigenvectors for the third member of the triplet have i replaced by —i, ie they are as if 
the complex conjugate was taken with the parameter(s) being regarded as real. 



Curiously the eigenvectors for the eighth triplet are precisely the same as for the third 
triplet. 

The importance lies in the fact that for each of the combinations of exponents of the 
leading-order terms and their resonances the number of arbitrary constants is the correct 
number. 



4 Interpretation 



The first triplet presents Right Painleve Series as the solution. In addition to the arbitrary 
location of the singularity three constants enter at the resonance +1 and two at the resonance 
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+2. The first few terms of the Laurent expansion for the first member of the triplet are 

X = |t"^ + ao - {61 + ci + (ao + 60 + co)^} T + I |(ao - ^o)^ + (co - ao)^ - 3ao(6o + Co) 

+ 3{bi + ci){ao-bo}T^ + ... 
y = Ir'^ + bo + biT-^ |coCi + Cgao - bocl - Cibo - 2aoboCo - blco - 2aobl - b^ + a^bi - b^bi^r'^ + ... 
z = ^T^^ + Co + ciT + I |cq + 2aoCo + 6oCo + 2ao&oCo + blco - ao&o + co(&i + ci) - aoCi - 60^1} + 

It is easily appreciated that the coefficients become unmanageably complicated within a 
few terms. Nevertheless it is evident that three arbitrary constants, (cq, bo, Cq), enter at 
the first resonance and two arbitrary constants, (61, Ci), at the second resonance. Moreover 
there is consistency at the second resonance and the arbitrary constants introduced at the 
first resonance remain arbitrary. Altogether the first triplet belongs to the so-called usual 
case of the singularity analysis with the single generic resonance at —1 and multiple positive 
resonances with the fortunate property that their geometric multiplicity equals their algebraic 
multiplicity at each resonance. 



All the other possible patterns for the resonances are of mixed sign and give rise to Laurent 
expansions in which the exponent goes to plus and minus infinity and so each represents an 
expansion on an annulus centred on the singularity [1] . This is to be expected when there is 
a resonance at —2. However, it is not so expected in the third triplet where one only has a 
double 'generic' —1 resonance. Unfortunately there does not appear to exist an algorithmic 
procedure which enables one to calculate the coefficients of the series by substituting them 
into system (12. ip . Nevertheless we can see how two arbitrary constants enter at the resonance 
— 1 by partially following the example of Latifi et al |T5] to construct a particular solution 
by putting to zero the constants of integration which enter at some of the resonances. As 
we are interested only in the descending series, we set to zero the three arbitrary constants 
occurring at the resonance +1 and one arbitrary constant entering at +2, ze, we write the 
solutions in the form of Left Painleve Series. We choose the first member of the triplet to 
illustrate our point. We obtain 

2 02 02 - 20262 - 2^2 9(a2 + 62)0262 
X = ^ + — + -3 1 + ••• 

1 ^ 62 ^ a2 + 40262 + bl ^ 3(4 + 30^62 - 6^) ^ 

3t 

1 02 + 62 02 + 0262 - 62 3(a^ - 30262 - 



3t 



+ 



We see that the expansion is in terms of two arbitrary constants, 02 and 62. It might 
be thought that there are three arbitrary constants since we also have the location of the 
movable singularity. However, if we expand the first term in each series in the usual way to 
explain the existence of the generic resonance, we see that to is absorbed into 02 and 62 
We emphasise that we took the particular values of the constants of integration made above 
so that we could provide an accurate representation of the solutions and the way in which 
the two arbitrary constants entered into the series at the double —1 resonance. The general 
solution must have the arbitrary constants from the positive resonances as well as the two 
entering at the negative resonance. 
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In f l3.2l) . which we recall is common to both the first member of the third triplet and the 
first member of the eighth triplet, the eigenvector corresponding to the triple —1 resonance 
is 




/U + 





(4.1) 



ie, there are three arbitrary constants entering at this resonance. As in the case of the first 
triplet the location of the movable singularity is absorbed into these three constants so that 
system ( 12. ip has six arbitrary constants with the other three entering at the double —2 and 
single +1 resonances. 



5 Conclusion 



We have seen by way of explicit example on the system, fl2.ip . of three nonlinear second- 
order ordinary differential equations that the so-called generic resonance at —1 plays a more 
complex role in the singularity analysis than it is generally accorded in the literature. The 
reasons for this to occur are unknown and perhaps it would be better not to speculate. 
Speaking of our own experiences our understanding of the singularity analysis has evolved 
in time from the examination of many explicit examples. In retrospect one wonders how 
one could have failed to appreciate that the essential feature of the singularity analysis is 
the determination of Laurent expansions as solutions to the differential equations under 
consideration. Then one may have recalled that the expansion of a function as a Laurent 
series can be at various parts of the complex plane. It can be in the vicinity of the singularity, 
which is where we have a Right Painleve Series. It could be beyond the last singularity and 
so have the nature of an asymptotic expansion which is where we have a Left Painleve Series. 
Finally it could be between two singularities separated from the singularity about which the 
expansion is being made. This gives a series in which the exponent can range from minus 
infinity to plus infinity. 

It is with the appreciation of the different domains of convergence of the results of one's 
analysis that one understands the three possible types of series. That one obtains a series 
differing in appearance |1] from what the 'traditional' analysis allows does not mean that 
one has, for example, an essential singularity. There may be an essential singularity, but the 
evidence provided by the standard singularity analysis is not going to enable one to state 
that as a certainty as was the case for the Mixmaster Universe [S]. 

The value of system (12. ip is that we have the solutions in closed form. Consequently we 
can look at the results of the analysis and see whether they say anything peculiar about an 
integrable system. Moreover, when we obtain similar results from a system the integrability 
of which is unknown, we must be careful not to make an inference which is unsupported by 
the evidence. 

Finally we have seen that a multiple resonance at —1 can indicate the entry of more 
than one constant into the solution. This is important in a consideration of integrability. 
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Perhaps an unexpected result was that a multiple resonance at —1 without any other negative 
resonances could presage Laurent series going to minus infinity. 

There is a question which remains unanswered. The Laurent series which arise as a result 
of the singularity analysis can be expected to have limitations on the regions in which they 
are convergent. If one has the three possible types of series, one would expect that the whole 
of the complex plane would be covered apart from the singularities. What happens, as we 
have in the case of (12.11) . when the series are convergent on the interior of a punctured disk 
and the interior of an annulu^ centred on a singularity? There are no asymptotic series (Left 
Painleve Series). Is it necessary to cover the whole complex plane, apart from singularities, 
to be able to infer integrability? 
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Coefficients {a, (3, 7} 


Resonances 


Triplet 


{0,0,0} 






(1 1 1\ 

13' 3 ' 3 j 

{i,i(-l + .x/3),i(-l-^x/3)} 


-1, 1(3), 2(2) 
-1, 1(3), 2(2) 
-1, 1(3), 2(2) 


1 


f2 2 2I 

l3' 3' 3/ 

|i,U-i-i^),ir-i+iv^n 

-l + i^),|(-l-iV^)} 


-2, -1, 1(2), 2(2) 
-2, -1, 1(2), 2(2) 
-2, -1, 1(2), 2(2) 


2 


(2 _i 

13' 3 ' 31 

{i,l(l-iV3),i 1 + .^)} 
{i,|(l + zV3),|(l-i^)} 


-1(2), 1(3), 2 
-1(2), 1(3), 2 
-1(2), 1(3), 2 


3 


{1,0,0} 






{l,i(-3-z^V|f-3 + 2y3)} 
{l,i(3-zx/3) ,i(3 + zy3)} 


-2, -1(2), 1(2), 2 
-2, -1(2), 1(2), 2 
-2, -1(2), 1(2), 2 


4 


{1, |(-3 + i^K U-3 - za/s) } 
{l,i(3 + ^^),i(3-^V3)} 


-2, -1(2), 1(2), 2 
-2, -1(2), 1(2), 2 
-2, -1(2), 1(2), 2 


5 


|4 _2 2| 

|l,Ui-'v3U(i+.v/3n 

{|,|(l + zV3),|(l-iV^)} 


-2(2), -1(2), 1, 2 
-2(2), -1(2), 1, 2 
-2(2), -1(2), 1, 2 


6 


-l + ^V3 ,i -l-^V3 } 


-2, -1(3), 1(2) 
-2, -1(3), 1(2) 
-2, -1(3), 1(2) 


7 


r 5 1 1 1 

|3'-3'-3} 
||,iri-ZV^),iri+2V^U 

{|,i(l + ^V3),i(l-^^/3)} 


-2(2), -1(3), 1 
-2(2), -1(3), 1 
-2(2), -1(3), 1 


8 


{2,0,0} 
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